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Abstract 

We study the N-dimensional pressureless Navier-Stokes-Poisson equations with density- 
^>f-^ ' dependent viscosity. With the extension of the blowup solutions for the Euler-Poisson equa- 

' tions, the analytical blowup solutions, in radial symmetry, in J?^ (A'^ > 2) are constructed. 



1 Introduction 

The evolution of a self-gravitating fluid can be formulated by the Navier-Stokes-Poisson equations 
of the following form: 



X 

: I Pt+v • (pu) =0, 

(pu)j+V • (pu (g) u) + VP ^-pV* + vis{p, u), (1) 
A$(t,x) = a{N)p, 

where a{N) is a constant related to the unit ball in : a{l) — 2; a(2) = 27r and For > 3, 

aiN) ^ NiN - 2)ViN) ^ NiN - 2) ^^^^^ ^ , (2) 

where V{N) is the volume of the unit ball in and F is a Gamma function. And as usual, 
p — p{t, x) and M = {[{t, x) G are the density, the velocity respectively. P = P{p) is the 
pressure. The evolution of the cosmology can be modelled by the dust distribution without pressure 
term. That describes the stellar systems of coUisionless, gravitational N-body systems. 
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In the above system, the self-gravitational potential field $ = x) is determined by the 
density p through the Poisson equation. 

And vis{p,u) is the viscosity function: 



vis{p,u) := v(m(p) V •")• 



(3) 



In this article, we seek the radial solutions 



X X 

p{t, x) = p{t, r) and u = —V{t, r) :— —V, 

r r 



(4) 



with r = (^J2iLi ■ particular, we can deduce the viscosity function ([3]) to be 



V 



d d 



d 



dxi dx2 dx 



N 



/a^|_ X2_ 

\ r ^ r ^ 



Xn 



V 



N 



fe=i 



dxk V r 



2 2 

r — xt 



fe=i 



r r 



= V 
= V 



Kp) ( V + Vr 



N -I N -I 

= Vtl{p)i V + Vr) + /i(p)V( V + Vr) 

r r 
d / N -I \ d N -1 



N ~l\ dr _ N -I dr 



dx, 



■V 



Vr 



(5) 
(6) 

(7) 

(8) 

(9) 
(10) 
(11) 



dr 



dxi dx, 



-f^rip) (^V + V,.]+p{p) 

r \ r 



^Lr {P) (^V + Vr]+p{p) 



N — \ \ Xi ^ ^ N — 1 Xi^^ Xi 

-V + -Vr + — 



"^-'K+^Vr + Vrr 



Vr, 



(12) 



(13) 
(14) 
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The equations ([3]) can be converted to the radial solutions to have 



p [ut + (m • V) u] + VP(p) + 

dui \ dP{p) 



fc=i 



du 



dxi 
dui 



p {u^)t+n,— + E^k^ 



dx 



dxk 



9$ (p) 

^p-^jQri^l,2,...N 

dP{p) 



dxi 



dxi 



d_ 

dt 



\ r / r oxi \ r / r oxt ^ r / 



f dP{p) 9$ {p) \ dr_ 



(15) 
(16) 

(17) 
(18) 



-Vt 



-v 



Xi 



Vt 



k^i ' 



\ dr ' 'dr j dxi 

~XiX]^ XiXj^ \ r / \ / \T Xi 



VVr 



V + -^Vr + [Pr{p) + (p)] 



Prip)+p^r (p) 



- [p {Vt + VVr) + Pr[p) + P^r {p)] ■ 

r 



(19) 
(20) 
(21) 



Therefore, the Euler-Poisson equations in radial symmetry can be written in the following form 



Pt + Vpr + pVr + 



N ~1 



-pV = 0, 



P (Vt + VVr) + Pr{p) + (p) = Pr {p) {^V + Vr) + p{p) [- (^) V + + K, 



Here we under a common assumption for: 



p{p) := Kp\ 



(22) 



(23) 



and K and 9 > Q are the constants. For the study of this kind of the above system, the readers 
may refer [7], [9], [10], [12] and [14]. In particular, when = 0, it returns the expression for the 
V-dependent only viscosity function: 



N -I 



Pt + Vpr + pVr + pV = 0, 

r 

p {Vt + VVr) + Pr{p) + P'^r {p) = p{p) [Vrr + " (^) V] 



(24) 



which are the usual form Navier-Stokes-Poisson equations. The equations ([T])i and ^2 {vis{p, u) ^ 
0) arc the compressible Navier-Stokes equations with forcing term. The equation ([Ija is the Poisson 
equation through which the gravitational potential is determined by the density distribution of the 
density itself. Thus, we call the system ([1]) the Navier-Stokes-Poisson equations. 

Here, if the vis{p,u) — 0, the system is called the Euler-Poisson equations. In this case, the 
equations can be viewed as a prefect gas model. For = 3, ([1]) is a classical (nonrelativistic) 
description of a galaxy, in astrophysics. See [2], [8] for a detail about the system. 

P = P{p) is the pressure. The 7-law can be applied on the pressure P{p), i.e. 



P{p) = Kp-< 



(25) 
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which is a commonly the hypothesis. The constant 7 — cp/cy > 1, where cp, Cy are the specific 
heats per unit mass under constant pressure and constant volume respectively, is the ratio of the 
specific heats, that is, the adiabatic exponent in (|^ . With K = 0, we call that the system is 
pressureless. 

For the 3-dimensional case, we are interested in the hydrostatic equilibrium specified by u = 0. 
According to ^2], the ratio between the core density p(0) and the mean density pfor6/5<7<2 
is given by 

^ =(—y{^)) (26) 



P(0) X ~ . 

where y is the solution of the Lane-Emden equation with 72=1/(7 — 1), 

viz) + -y{z) + y{zr = 0, yiO) = a > 0, y(0) = 0, n = (27) 
z 7 — 1 

and zq is the first zero of y{zQ) = 0. We can solve the Lane-Emden equation analytically for 



yanal{z) ^ < 



1 - n = 0; 

sin z 

, n = 1: 



' 1 

= , 71 = 5, 



(28) 



and for the other values, only numerical values can be obtained. It can be shown that for n < 5, 
the radius of polytropic models is finite; for tt, > 5, the radius is infinite. 

For the isothermal case (7 = 1), the corresponding stationary solution is the Liouville equation, 

2 a(3) ^'""^ 

y{z) + -y{z)+^e =0. (29) 
z K 

Gambin [4 and Bezard P] obtained the existence results about the explicitly stationary solution 
(u — 0) for 7 = 6/5 in Euler-Poisson equations: 



V 2n 
where A is constant. 

In the following, we always seek solutions in radial symmetry. Thus, the Poisson equation 
is transformed to 

r^-i$„ {t,x) + {N -l)r^-Hr^a{N)pr^-\ (31) 

^r = ^j^ r Pit,s)s''~'ds. (32) 
Jo 

In this paper, we concern about blowup solutions for the A^-dimensional pressureless Navier-Stokes- 
Poisson equations with the density-dependent viscosity. And our aim is to construct a family of 
such blowup solutions. 

Historically in astrophysics, Goldreich and Weber [5] constructed the analytical blowup solu- 
tions (collapsing) of the 3-dimensional Euler-Poisson equations for 7 = 4/3 for the non-rotating 
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gas spheres. After that, Makino [B] obtained the rigorously mathematical proof of the existence 
of such kind of blowup solutions. And in [5], Deng, Xiang and Yang extended the above blowup 
solutions in {N > 4). Recently, Yuen obtained the blowup solutions in with 7 = 1 in [13] . 
The family of the analytical solutions are rewritten as 
For TV > 3 and 7 (2iV - 2)/N, in [13] 

1 



pit, r) = < 



ait) 



WViW)^'"^''-'^^ for r<a(t)Z 



0, for a{t)Z^, < r. 



N ci(t) 
a{t) 



(33) 



y(^)+ yiz)+- 



:y(z)^/(^-')= M, y(0) - a > 0, 2)(0) = 0, 



z - • {2N -2)K- 
where /i = [N{N - 2)A]/(27V - 2)ii: and the finite Z,, is the first zero of y{z); 
For TV = 2 and 7 = 1, in [l4| 

1 e^('-M*)),F(t,r)=^r; 



A 



ait) 



a{t)= -— , a(0) = > 0, a(0)= ai; 

a(2) 



(34) 



2/(^) + -2/(2)+ 



Ai, y(0) = a, y{0) = 0, 



where if > 0, = 2X/K with a sufficiently small A and a are constants. 

For the construction of special analytical solutions to the Navier-Stokes equations in i?^ without 
force term, readers may refer Yuen's recent results in |15j . 

In this article, the analytical blowup solutions are constructed in the pressureless Euler-Poisson 
equations with density-dependent viscosity in R^ in radial symmetry: 



iV - 1 

Pt + Vpr + pVr + pV = 0, 

r 

piVt+VVr) 



a{N)p 



r.N-1 



p{t, S)s^-^ds + [KpO]r {^V + Vr) + {Kp^)iVrr + 



TV- 1, 



(35) 



in the form of the following theorem. 



Theorem 1 For the N -dimensional pressureless Navier- Stokes- Poisson equations in radial sym- 
metry, ^35\). there exists a family of solutions: 
for N>2 and N with 9 = (2TV - 3)/TV, we have, 

1 



Pit,r) 



(T - Ct)i 



Ny\T-Ct 



for 



T-Ct 



< Zn 



, V{t,r) 



-C 



0, for Zo<^ 



T-Ct ' 



T-ct ■ 



,s N^l a(N)(N-'S) 



z - ' ■ (2TV-3)CTVk 
for TV = 3, with 0=1, we have, 

1 



< 



Pit,r) = 



-C 



T-Ct 



[T - Ct)3 
2 47,- vi^) 

y{z) + -y{z)+-—e = 0, 2/(0) = «, m = 0, 

Z O i V 



(36) 



(37) 
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where T>0, k>0, C^O and a are constants, and the finite Zq is the first zero of y{z). In 
particular, for C > 0, the solutions blow up in the finite time T/C. 

2 Separable Blowup Solutions 

Before presenting the proof of Theorem[l] we prepare some lemmas. First, we obtain the solutions 
for the continuity equation of mass in radial symmetry psp i . 

Lemma 2 For the N- dimensional conservation of mass in radial symmetry 

N -1 

Pt + Vpr + pVr + pV = 0, (38) 

r 

there exist solutions, 

, ^ fir/iT-Ct)) , , ~C 

(r-ct)^ (39) 

where f > E , T and C are positive constants. 

Proof. We just plug ^ into Then 

TV - 1 

Pt + Vpr + pVr + pV (40) 

r 



(-iV)(-C)/(r/(T ~ Ct)) hC)rf{r/{T - Ct)) 
(T - Ct)^+i (T - Ct)^+^ 

{-C)rf{rl{T^Ct)) , f{r/{T-Ct)) (-C) , N - I f{rl{T - Ct)) {~C) 



(41) 



(42) 



T-Ct {T~Ct)^+^ (T-Ct)^ T-Ct r {T - Ct)^ T - Ct 
= 0. (43) 

The proof is completed. ■ 

Besides, we need the lemma for stating the property of the function y{z) of the analytical 
solutions (|36p and ([57]) . We need the lemma for stating the property of the function y{z). In 
particular, the solutions (j36[) in 2-dimcnsional case involve the following lemma. The similar lemma 
was already given in Lemmas 9 and 10, in |14j . by the fixed point theorem. For the completeness 
of understanding the whole article, the proof is also presented here. 

Lemma 3 For the ordinary differential equation, 



1 

—1 

z' 

y{0) = a > 0, 2/(0) = 0, 



m+-yi^) - ^ - ^^^^ 



where a is a positive constant, 
has a solution y{z) € C^ and lim y{z) = oo. 



2 — * + 00 
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Proof. By integrating (HU, we have, 

u rz ^ 

y(^) = - / -T^-""^-' ^ 0- (45) 
Thus, for < z < zq, y{z) has a uniform lower upper bound 

V{z) > 2/(0) = a > 0. (46) 

As we obtained the local existence in Lemma [31 there are two possibilities: 

{l)y{z) only exists in some finite interval [0, zq]: (la) lim y{z) = oo; {lh)y{z) has an uniformly 
upper bound, i.e. y{z) < ao for some constant ao. 

{2)y{z) exists in [0, +cxd): (2a) lim y(z) = oo; (2b)i/(z) has an uniformly upper bound, i.e. 
y{z) < P for some positive constant /3. 

We claim that possibility (1) does not exist. We need to reject (lb) first: If the statement (lb) is 
true, (|45l) becomes 

Thus, y{z) is bounded in [0, zq]. Therefore, we can use the fixed point theorem again to obtain a 
large domain of existence, such that [0, + for some positive number 5. There is a contradiction. 
Therefore, (lb) is rejected. 

Next, we do not accept (la) because of the following reason: It is impossible that lim y{z) = oo, 

z^zo- 

as from (1471) . y{z) has an upper bound in [0, zq]: 

S > yi^)- (48) 

Thus, (gB]) becomes, 

"° <yzo a{zo) 



y{zo) = y(0) + / yis)ds <a+ j ^ds = a + (49) 



Since y{z) is bounded above in [0, zq], it contracts the statement (la), such that lim y(z) = oo. 

z-*zo- 

So, we can exclude the possibility (1). 

We claim that the possibility (2b) doesn't exist. It is because 



Then, we have. 



Viz) ^ds>^- ^ds = ^. (50) 
z Jo V{s) z Jq (}^ 2f3^ 



yiz)>a+—z\ (51) 



By letting z — > oo, (|5ip turns out to be, 

y{z) = oo. (52) 

Since a contradiction is established, we exclude the possibility (2b). Thus, the equation (|^^ exists 
in [0, +oo) and lim y{z) = oo. This completes the proof. ■ 

z — ^-\-oo 
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For N > A, our blowup solutions depend on the Lane-Emden equation, 



iV - 1 . 

y{z) H — —yi^) + ^y(^)'= 0' 



(53) 



y{0) = a > 0, m = 0, 

where s > 1 and a are positive constants, which is a particular case of the Emden-Fowler equation 



il + z^-"u" = 0, 



(54) 



with n > 1, after the transformation 



(55) 



For the existence and uniqueness of the equation (j53p . the interested readers may refer the survey 
paper [11]. 

The proof of Theorem [1] is similar to the ones in [3], [13] and [Mj. The main idea is to put the 
analytical solutions to check that if they satisfy the system ((35]) only. 
Proof of Theorem [Jl From Lemma [21 it is clear for that (|36| satisfy (|35]) i . 
For the case of iV > 2 and N ^ 3 with 6 = {2N - 3)/N, we plug the solutions ^ into the 
momentum equation (j35[) 9. 



(_C)(-1)(-C) , i-C) i-C) 



or 



[T - Ctf 
1 



[T-Ct)^" \T-Ct 



T-Ct T-Ct 

N/{N-3)-\ (■^N-3)/N 



' s as 



r"~^ J [T-Ct) 





AT 



T^Ct 



(56) 
(57) 

(58) 



rN-1 J {T- Ct) 





N 



\T~Ct ) 



TV - 3 {T-Ct) 



N-3 



(T - Ct)^+^ 



c^{N)p /-ylr^)^/^^-'^ ^ {2N - 3)CNk p 



r^-i J {T~Ct) 





N 



N-3 {T -Ct)^-^" \T -Ct 



(T - Ct) 



N-l 



N-3 y^T-Ct^^ rN-^{T~Ct)J^^T~Ct' 





(T-Ct)^-i 
P 



(2N - 3)CNk ■ , r , a{N) 



vliT-Ct) 



N-3 '''T-Ct' (=^)JV-i 



y{s)^/iN~3)s^-^ds 



-Q 



(T-Ct)^-^^ \T-Ct 



(59) 
(60) 

(61) 

(62) 
(63) 
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And denote 





Differentiate Q{z) with respect to z, 



(64) 



?/(z)+Q;(7V)y(z) 



N -3 



y{s) 



N/{N-3)N^1 



s^'-'ds (65) 



iV-3 



-(N-l) 



(66) 



(67) 



where the above result is due to the fact that we choose the following ordinary differential equation, 

(68) 



■■/ ^ a{N){N-3) , .N/tN „ 



z ^ ■ ■ (27V-3)C7Vk' 
y(0) = a>0, y(0)=0, 



in Lemma|3]and the well-known results in the Emden-Folwer equation. With Q(0) — 0, this implies 
that Q{z) — 0. Thus, the momentum equation (|35p 9 is satisfied. 

For the case of iV = 3 with 9 ~ 1, we plug the solutions ([37| into the momentum equation ([35]) 2, 



N - 1 



piVt + VVr) + / pit, S)s'ds - [p{p)]r V +Vr]- fiipWrr + -Vr " -^V) (69) 



(-C)(-l)(-C)^ ^ (-C) ^ (-C) 



(T - Ci)2 



47r/3 /" e' 



,y(s/(T-Ct)) 



r2 7 (T - Ct)3 



-s ds 



/ N , , / 2 (-1) 2 (-1) 



„y(r/{T-Ct)) 



1 NC Anp /•e2'(^/(^-c^*)) 



(T-Ct)3 "'T-Ct' T -CtT -Ct r"^ J [T - Ct) 



-s ds 



CNkp ., r ^ Anp z" e^'^''/^^-^*)) 



• yirr /-r. ) „2 



[T -Cty^T ~Ct' r2 y (T-Ct) 



-s (is 



(T - C<)2 



CNny{——) + 



47r 



T-C<' r^{T-Ct) 



(70) 
(71) 
(72) 

(73) 

(74) 



By letting uj = s / {T - Ct), it follows: 



[T - Ctf 



CNKiji——) + 



r/{T-Ct) 



T-Ct (TTct)^ 



(75) 



10 



M.W.Yuen 



And denote z = r/{T -Ct), 



'T-Ct 

Differentiate Q{z) with respect to z, 



Q{z) = CiV«;y(z)+47re^(^) + ^ / e^^'^) 



2CN 



z 

wliere tfie above result is due to tlie fact that we choose the Liouville equation: 

2 47r 



(77) 

(78) 

(79) 
(80) 

(81) 



y(0) = a, y(0)= 0. 

With (5(0) = 0, this implies that Q{z) = 0. Now we are able to show that the family of the 
solutions blows up, in the finite time T/C. This completes the proof. ■ 
The statement about the blowup rate will be immediately followed: 



(82) 



Corollary 4 The blowup rate of the solution t36\) and |j'7| ) is 

lim p(t,0)(T- C<)^ > 0(1). 



T_ - 

C 



Remark 5 Besides, if we consider the Navier-Stokes equations with the repulsive force in radial 
symmetry in (N 1), 



iV - 1 

Pt + Vpr + pVr + pV = 0, 

r 

p{Vt+VVr) 



a{N)p 



rN- 



f / p{t, S)s^-^ds + [kP% {^V + Vr) + i^p'Wrr + K 



iV-1,. A^-1 



V) 



(83) 



the special solutions are: 

for N>2 and N ^3, with 9 = {2N - 3)/iV, 

1 .JL. 

Pit, r) ^ { i^^^ClW'^^'^ ^-^,for^,<Zo; 



V{t,r) = 



-C 



0, for Za < 

T— I 



T-Ct 



ct 



(84) 



A^- 1 ^ a{N){N -i) ^ 



y{z)+^y{z) (2A.-3)CA^^c 
the finite Zo is the first zero of y{z); 
for iV = 3, with B=l, 

1 



y{z) 



P{t,r) 



0, y(0) = a > 0, 2/(0) = 0, 



C 



^y{r/{T-Ct))^ ^ 



T-Ct 



[T - Ct)^ 

2 An y^'^^ 

y(z)+-y{z) ~ = 0, y{0) = a > 0, y(0) = 0. 



(85) 
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